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We present a realistic supergravity inflation model which is free from the overproduction of poten- 
tially dangerous relics in cosmology, namely moduli and gravitinos which can lead to the inconsis- 
tencies with the predictions of baryon asymmetry and nucleosynthesis. The radiative correction 
turns out to play a crucial role in our analysis which raises the mass of supersymmetry breaking 
field to intermediate scale. We pay a particular attention to the non-thermal production of grav- 
itinos using the non-minimal Kahler potential we obtained from loop correction. This non-thermal 
gravitino production however is diminished because of the relatively small scale of inflaton mass 
and small amplitudes of hidden sector fields. 
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I. INTRODUCTION 

There exist generic and well known problems in con- 
structing supergravity inflation models with broken lo- 
cal supersymmetry in vacuum. First, one must carefully 
choose the superpotential and Kahler potential such that 
non-renormalizable terms do not spoil the flatness of the 
inflaton potential, which is also essential to obtain the ob- 
served CMB spectrum. Second, most supergravity infla- 
tion models run into cosmological problems at late times 
due to an over-abundance of harmful relics such as mod- 
uli and gravitinos, and we need to check that a model is 
free from these problems to be consistent with the data 
of baryon asymmetry and nucleosynthesis predictions. 

Several supergravity inflation models free from abun- 
dant moduli have been proposed [1-5], and we give a brief 
review for this moduli problem with emphasis on the im- 
portance of the radiative corrections. Especially when 
the supersymmetry breaking field has flat direction at 
tree level, radiative correction has a significant effect on 
its potential to lead to the necessity for the modification 
of minimal Kahler potential. 

More recently, the nonthermal production of graviti- 
nos has been drawn attention [6-10]. These problems 
on gravitino production during preheating era, however, 
have been analyzed only using the minimal Kahler po- 
tential so far. The gravitino interactions depend on 
the form of Kahler potential and consequently the non- 



thermal production of gravitinos can depend on its form 
as well. Our treatment is the first analysis of grav- 
itino non-thermal production taking account of the non- 
minimal Kahler potential obtained from the loop correc- 
tion of supersymmetry breaking field. 

The detailed analysis for nonthermal production of 
gravitinos in a system of coupled fields has been done 
only for models with non-renormalizable hidden sector 
supersymmetry breaking, and only the Polonyi model in 
particular [7,11]. It was shown, in this model, that dom- 
inant fermion fields which arc created efficiently through 
preheating mechanism are inflatinos rather than grav- 
itinos (thus free from gravitino problem). This type of 
supergravity inflation model is a good toy model to inves- 
tigate the nonthermal particle production in preheating 
era, but is not a realistic inflation model in that it sill suf- 
fers from Polonyi problem which is as serious a problem 
as gravitino problem. Wc give the first realistic super- 
gravity inflation model in this sense where nonthermal 
production of gravitinos as well as moduli problem were 
explicitly analyzed. 

The paper is structured as follows. In section II, we 
explain our choice of superpotential and discuss how in- 
flation develops in our model. We then discuss how su- 
persymmetry breaking field evolves and calculate its ra- 
diative correction and its modification to minimal Kahler 
potential arising from this loop correction. In section 
III, we see if our model leads to any cosmological cri- 
sis, namely, moduli and gravitino problems. We give the 
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conclusion and discussion at the end. 



II. SETUP 

We consider the superpotential [1] consisting of the 
inflaton sector and hidden sector of O'Raifeartaigh type 
[13] which are gravitationally coupled to each other. 



W 



(E - Mf 



+ $1 (k$2 - At 



A$2*3 + C, 

(2.1) 



where supcrficid S includes inflaton scalar component 
a and $i includes O'Raifeartaigh scalar field (/>i. A ~ 
10~*Mp from COBE normalization and M is set to 
Mp{= 2.436 • lO^^GcV) so that the inflaton potential 
keeps the flatness around the origin (i.e. for -^V{0) ~ 

0, -^V^O) ~ ). The dimensionless parameter k is of 
order unity while the other mass parameters fi and A are 
of intermediate scale ( ~ 10~^ Mp ). C is the constant 
term to cancel the cosmological constant at the vacuum. 
We start with the discussion for the evolution of scalar 
fields in infiaton sector and hidden sector which are only 
gravitationally coupled to each other and, for the sake of 
clarity, wc first treat each sector separately followed by 
the discussion including the coupling with non-minimal 
Kahler potential arising from the radiative correction. 
We assume gauge singlets in the potential for simplicity 
in the following. 



potential from Winfuton for the real part of scalar com- 
ponent, a, becomes ( in natural units ) 




(2.6) 



(2.7) 



We here point out the absence of linear and quadratic 
terms which enables the potential to keep the fiatness 
around the origin. We can obtain the value of A ~ 
10"'* Mp from COBE normalization condition [15] *, 



~ 0.027Mp(l - 3.2e + 0.5r]), 



(2.8) 



which should be evaluated at the horizon exit. Scale of 
inflaton flcld when the inflation ends and the cosmolog- 
ical scales leave the horizon are obtained from slow-roll 
conditions, e < 1, ry < 1, and 60 e-folding condition. 



60. 



(2.9) 



We also note the scale of inflation is of the order A^ and 
the mass of the inflaton is of the order A^/Mp with its 
decay width V„ ^ m^^/Mp = A^/Mp assuming gravita- 
tional strength coupling to ordinary fields. 



B. Hidden Sector 



A. Inflaton Sector 

The superpotential of inflaton sector is given as 
.2 i^-Mpf 



Winflaton = A 



(2.2) 



p 



The general expression for supergravity potential be- 
comes, for Kahler potential, K, and superpotential, W, 



with 



V = mi K-^/ m? - 3 M-Vl^ 



K'j= „ , „ , , m = e^W 



,1 _ T-lI 



D'm = d'r 



1 



2M2 



(2.3) 



(2.4) 

(2.5) 



If we consider, for the moment, the minimal form of 
Kahler potential, K = T,^ Y,, the effective supergravity 



The supcrsymmetry breaking sector is that of 
O'Raifeartaigh model, 



idden 



$1 (K$i - /i^) 



A$2$3 + C 



(2.10) 



This is a familiar example of supcrsymmetry breaking 
due to non- vanishing F-term from $i, \F\ = /j."^, in the 

2 

vacuum. Therefore we add C = ^Mp (compare with 

term in eqn(2.7) ) for the vanishing cosmological 
constant at the vacuum.'^ We should, however, expect 



* [1] gives an order of estimates lO^^Mp < A < W'-^ -'Mp 
from the constraints on gravitino abundance and proton 
decay. 

^This additional constant term C, strictly speaking, should 
be modified if we include the radiative correction and the 
coupling between inflaton and hidden sectors. We, however, 
stick to this value of C for simplicity because this modification 
essentially does not change our discussion. 
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that, when the fields are far away from the vacuum dur- 
ing the inflation, there are additional .F-terms from other 
fields in the effective scalar potential and these i^-terms 
lead to additional 'cosmological constant' A* [3]. Adding 
e^A^ in the potential indicates us that this cosmologi- 
cal constant term during the inflation gives an additional 
effective mass of order to each field in the model. 
Fields (f>2 and cf)^ do not posses the linear terms and these 
two fields quickly roll down to the origin ( i.e. to their 
minimum) during the inflation. The scalar field , how- 
ever, has a liner term and its minimum shifts according 
to the evolution of infiaton field as we shall see in the 
next section. Because we are interested in the particle 
production after the infiation, we focus on the evolution 
of O'Raifeartaigh field 0i among the fields in this super- 
symmetry breaking sector. Moreover, the F-term at the 
vacuum has the contribution only from the scalar field 
which turns out to have fiat direction at the tree level. 
Because of this flat potential with respect to at the 
tree level, the global supersymmetry radiative corrections 
have a significant effect on the effective potential and 
consequently give non-negligible modification to the min- 
imal Kahler potential. The radiative corrections of local 
supersymmetry are always Planck mass suppressed and 
we do not consider them here. We calculated this non- 
minimal Kahler potential from the calculation of loop 
correction [16,17]. Setting the parameter range to be 
2k/x^ < to make 02 and (ps stay at the origin in the 
vacuum leads to the following one-loop correction due to 



(2.13) 



one loop 
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(2.11) 



(2.12) 



where we have defined 

Mf = 1(^1 - A,), M| = l(Ai + A^), 

Mi - ^(^3 - ^4), M| = ^(^3 + A^), 

Nl = \{B^ - B^),Ni = \{B^+ B^), 

^1 = 2A^ - 2K/i2 + 4K2|0l|^ ^3 = 2A^ + 2Kfi^ + 4:K^\(1)i\ 

A2 = ^^K^^l'^ + 16A2/v2|,/,i|2 - 16k3^2|,^i|2 + 16Av4|0i|4, 

A^ = ^AK^fi^ + 16A2k2|(/)i|2 + 16kVI0iP + 16k4|0i|4, 
Bi=2X^ + B2 = v'16A2k2|^i|2 + 16k4|0i|4 . 



We have used the MS scheme and taken the renormal- 
ization scale to be A. We are concerned with the regime 
101 1 < A/k, as we shall show in the next section. In this 
case, and for 2kij,'^ <C A2, we can approximate the one 
loop potential as 



Vr. 



one loop 
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A2 



(2.14) 



where . . . are terms of order k** 1 0i | ^ / A^ and higher and C\ 
is a small constant term (^ /x'*) which can be absorbed 
into the constant part of the superpotential. We can find 
the following one loop correction to the Kahler potential 
{K = K^.nimai + Kcorrection), by Comparing eqn(2.14) 
with eqn(2.3), 



K. 



correction 
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A2 



(2.15) 



We shall use this non-minimal Kahler potential in our 
analysis for non-thermal production of gravitinos. We 
note that this radiative correction enhances the coupling 
to the longitudinal component of gravitino and raises the 
mass of 01 which was massless at the tree level to the 
intermediate scale 

4 4 2 4 9 

"^01 = = ^ ^ith = (2-16) 

which turns out to be crucial to evade the moduli prob- 
lem. 



III. HARMFUL RELICS 

We are now in a position to discuss the fate of in- 
fiaton and O'Raifeartaigh fields in the coupled effective 
potential with non-minimal Kahler potential to see if our 
model leads to any cosmological crisis. We first briefiy re- 
view the resolutions of so-called Polon}d or moduli prob- 
lem and we further discuss the non-thermal production 
of gravitinos. 



A. Moduli Problem 

We here start with the discussion on well-known po- 
tentially dangerous problems, Polonyi problem or, in gen- 
eral, moduli problem. There are two aspects which we 
should worry about before our discussion on decay of 
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moduli into gravitinos. One is the case when the mod- 
uli decay very late ( i.e. during or after the nucleosyn- 
thesis) which can jeopardize nucleosynthesis predictions 
because of ultra-relativistic decay products directly from 
moduli fields destroying the light elements, in particu- 
lar, ^He and D nuclei. The other is when the entropy 
release due to its decay is so big that it can over-dilute 
the baryon asymmetry well below its acceptable amounts 
(so-called 'entropy crisis'). Our model does not have ei- 
ther of these problems because the radiative correction 
raises its mass to as much as intermediate scale. Its de- 
cay width is indeed enhanced up to F^^ ~ m\ /\F\'^ ~ 



a^/Lt with = K;^/47r ~ 0(10 ^) and this is 



of orderlO~^^Mp. So O'Raifeartaigh field decays around 
IQ^^Mp^ in our model which is much before the nucle- 
osynthesis starts around ~ 10^° and even well be- 
fore the reheating starts due to the inflaton decay around 
1/r,, ~ MpVA^ ~ lO^^Mp-i. We however need a great 
care about the possibility of decay products with long 
life-time, especially gravitinos. Gravitino decay rate is 



of order r„ 



lO-^^Mp and 



its relativistic decay products, especially ultra-relativistic 
photon/photino, can destroy the light elements in nucle- 
osynthesis (photo-dissociation process) as we just men- 
tioned. The possible abundant gravitino production from 
O'Raifeartaigh fields can be caused by the energy release 
stored during the inflation by the shift of the minimum of 
SUSY breaking field potential as infiaton evolves. If this 
energy release is too big, it could lead to large amount 
of gravitinos and upset the nucleosynthesis predictions. 
We can see this is not the case for our model as fol- 
lows [1,4,20]. During the inflation, due to the couphng 
to the inflaton (cr ~ Mp), O'Raifeartaigh field amplitude 

2 

is around the intermediate scale of order ~ 
at the minimum of its potential. Therefore we can esti- 
mate the energy stored in this O'Raifeartaigh field when 
it starts oscillation ( i.e. t,j,^ ~ "^^i^) ^® most of 
order 



TO. 



(3.1) 



and its number density n,p^ in this coherently oscillating 
O'Raifeartaigh field is at most 



(*0l) 



(3.2) 



We can now estimate its number density to entropy ratio 
at the time of reheating for the gravitinos through the 
decay of {aXt = tr, say ) in an adiabatically expand- 
ing universe. Assuming, for the upper bound, solely 



decays into gravitinos with 100% branching ratio and us- 
ing s ^ ^g- g^r""* (with effective degree of freedom) and 
a? ~ for matter domination [21] due to the coherently 
oscillating inflaton field which dominates the energy in 
the universe, 



«3/2 



S{tr) 



0.445,r3 



RH 



This can lead to the estimate of after reheating 

by substituting 1.66 • T^jj^/Mp ^ H ^ f-^ for U, 



0.44M2 



O.44to0,A4 



We can now compare this value with one corresponding 
to the gravitinos produced by the scattering in the ther- 
mal bath in reheating era obtained in MSSM [18,19], 

n/n™d(T « IMeV) 1.1 • 10"" (j^^) • (3-5) 

Using s = 1.8 • g^firad and (<C MeV) ~ 3.36, we obtain 



n/s ~ 1.8-10- 



/ Trh 
\ lOiOQeV 



(3.6) 



flow of the gauge coupling. We can now transform 
eqn(3.4) by substituting (2.16) to the following form. 

This is smaller than the thermal production of grav- 
itino (3.6) by two orders of magnitude. The radia- 
tive correction therefore induces intermediate mass scale 
for O'Raifeartaigh field and it consequently makes the 
O'Raifeartaigh field energy released through the decay 
into gravitinos after the inflation small enough to evade 
the abundant gravitinos. Hence our model does not suf- 
fer from moduli problem as far as the constraint from 
thermal gravitino production is satisfied. 

We mention that non-adiabatic production of moduli 
fields in pre-heating era could lead to abundant graviti- 
nos [25] . We are, however, not concerned about the para- 
metric resonance effects for moduli fields because scalar 
coupling terms in Lagrangian in our model are trilinear 
in hidden sector fields, and those couplings to inflaton 
field have always Planck mass suppression [5]. 

The exception where this Planck mass suppression 
does not occur and preheating effects could be impor- 
tant is the coupling involving longitudinal components of 
gravitino, which is the subject in the following section. 
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B. Non-thermal Production of Gravitino 

It has been argued recently that parametric resonance 
mechanism in preheating era for the creation of gravitino 
can be much more efficient than the thermal one [6-10]. 
In this nonperturbative mechanism, the gravitinos can 
be created non-adiabatically through the amplification 
of vacuum fluctuation via rapid energy transfer from co- 
herently oscillating inflaton field which still dominates 
the energy density in the universe just after inflation and 
before the reheating era. 

In analyzing the gravitino field equations in the fol- 
lowing, we see that the equations for transverse and 
longitudinal components of gravitino decouple. While 
transverse component equation has a Planck mass sup- 
pressed coupling and thus gravitationally suppressed par- 
ticle creation, longitudinal componc;nt equation is free 
from Planck mass suppression and it could lead to the 
abundant gravitino production well above the constraint 
from thermal production of gravitino. Indeed, gravitino- 
goldstino equivalence theorem states that the equation 
for gravitino longitudinal component can be reduced to 
the cqiiation of goldstino in global supcrsynimctry in the 
limit of weak gravitational coupling. This warns us that 
gravitino longitudinal components could lead to its effi- 
cient copious production without Planck mass suppres- 
sion. Our model however has a Planck scale amplitude 
for inflaton field after inflation, and it is not obvious if 
this naive intuitive picture analogous to the goldstinos 
in global SUSY applies here. Therefore we apply here 
the formalism developed in [6,7] to calculate the num- 
ber density of gravitinos created through the nonthermal 
process. 

We first need to describe the evolution of scalar fields 
and fermion fields and their interactions. It is convenient 
to work with, among other possible choices, the following 
rescaled quantities in our numerical analysis, 



A 



'3 . a ^ II 



Mr, 



X- ^ i-f H-H^P V- ^ 



" A' 

(3.8) 



where H is Hubble constant, V" is a scalar potential from 
eqn (2.3) with non-minimal Kahler potential obtained in 
(2.15), 



K = SSt + + $2$2 + *3*3 



327r2 



A2 



(3.9) 



In terms of these rescaled quantities, the equations of 
motion for coherently oscillating scalar fields </>(= a,^i) 
read 



dt'^ dt dS 



(3.10) 



We omit ~ in the following discussion as long as it is 
clear from the contexts. We can concentrate on the field 
equations for a and (t>i because the other fields in super- 
symmetry breaking sector quickly roll down to the origin 
during the inflaton and stay there-'-. So we can let the 
amplitudes of 02 and (j)^ vanish after obtaining the equa- 
tion of motion for a and to see the field evolutions 
after the inflation. 

The Fermion equation follows from the supergravity 
Lagrangian 

e-^L = -^M^R- {d^ cj)') (9%) - V 
-^M^i^^Rf' + ^m i^^R T^-- V-.ii 

+ \m* 1'^" ^uL - [xj W + X' PXj] 

-m'^ Xi Xj - X' 
+ (2 /ST/V^^fl Y^' X^ d^<p, + ■ -iVL + h.c.) 
+ ( four fermion and gauge interaction terms ). (3.11) 

This Lagrangian includes chiral complex multiplets 
{<f)i, Xi) a^nd the Ricci scalar R. Subscript L and R denote 
its projection through operators Pl = {1 + 75)/2, Pr = 
(1 — 75)/2. Gravitino kinetic term shows up in the form 
of 



(3.12) 



R^ = e-\^''P^j5l.Dpij„, 
with covariant derivative 

£>M^. = ( + J<"7mn) St - . (3.13) 

The kinetic term for chiral fermion is 

D^Xi = (d^ + I'^r^mn^ Xi 

[djKd^^ -d^Kd^ci>j]xi + ri''xjd^<Pk (3.i4) 



4M2 



*Once these fields roll down to the origin, they stay at the 
origin to any higher order because of R-symmetry. 



5 



with Kahler connection Fl'' = K ^i'd^Ki^ and ^mn 
[7m,7n]/2 . Its mass term reads 



2Mf 

(3.15) 

The combination of matter fields gives left-handed com- 
ponent of goldstino 

VL = m'xi + {mx' K/. (3.16) 

The supersymmetry transformation of goldstino [6] 

5^ = -^(m^/, + H') e, ^ (3.17) 

indicates that gravitino mass and Hubble parameter sig- 
nal supersymmetry breaking. We can obtain the grav- 
itino equation from this Lagrangian, 

+ mip^ = (d^ - ^7^) 7^^,. (3.18) 

In solving this gravitino equation of motion, wo gauge 
away the goldstino ( unitary gauge ) and use plane- wave 
ansatz for the spatial dependence of f/ip 



Ak-x 



More- 
over it is convenient to decompose the space component 
of gravitino field into the transverse part '^i^ and trace 
parts 6 = 7*Vi and fc^Vi as 

= + (P^), e + {Pk)i kii^i, (3.19) 



where 



(P,), =-4-(3fc,-7i(fc,-7^)). (3.20) 

This leads to the the following succinct form of dy- 
namical field equations which describe the degree of free- 
dom corresponding to transverse and longitudinal com- 
ponents, 



7°ao-|-i7' h + 



0,7 



am 
4 



0, 



(3.21) 



(do + B + i^'ki'^'^ A\e fc2T = 0, (3.22) 

v J aa 

where 



m = PRm + PLm*, 



A = — (ai — 7° 02) 

a 



ai = 



= -K ( 



3H^ + 2H 



a = SMU H"^ 
3 



j) -^|mp, a2 = 2mt. (3.23) 



We can easily sec, reducing the equation into this form, 
eqn(3.21) describing the transverse component of grav- 
itino tjjf is decoupled from the longitudinal gravitino 
component, and its coupling to the other fields are Planck 
mass suppressed. So we hereafter pay our attention to 
the equation which describes the longitudinal component 
of gravitino, cqn(3.22). The form of T in (3.23) tells us 
that, in the absence of Kahler terms which mix the vari- 
ous left chiral superflelds, we need only worry about the 
fermionic partners of dynamical scalar fields. Further- 
more, for our supcrpotcntial, there is no mixing between 
the fermion associated with 0i and those of 02 and (j)^, 
as long as (f>2 = (1)3 = which is true because they stay 
at the origin due to R-symmctry once they roll down to 
the origin during the inflation. Thus, even though the 
effective masses of the fermions corresponding to 02 and 
03 are changing, those fermions do not contribute to the 
goldstino and we can concentrate on the evolution of the 
other fields for the purpose of our calculation. 

Based on the form of equation of motion involving 
two chiral superfields, we can infer the terms in the La- 
grangian which describe the interactions of the two fields 
of our interests, namely, 9 ( longitudinal component of 
gravitino) and T ( combination of chiral fermions orthog- 
onal to goldstino v). Those interaction terms lead to the 
equation of motion in the following matrix form, 



{j''do + i^'kiN + M)X = 0, 



(3.24) 



with vector X 
ized fields 



^i:^ consisting of canonically normal- 



2iy k. 



{a a 



0, 



2- ■ ^■ 



and diagonal mass matrix M is given by 



(3.25) 
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M = diai 

3 ma 



ma 
2MI 



2 VM2 



ai + aa2 



ma 
"2M2 



+ -^j^ai + aa2 + a{mu + m22) j (3.26) 



and matrix N 



-ai 



+ 1" 



-d2 -A 
-A d2 



(3.27) 



for at = Ui/a and A 



Q!o. The existence 



of off-diagonal terms warns us the non-trivial mixing of 
fermion eigenstates. 

Once we can reduce the Lagrangian into this form of 
matrix expression, we can obtain the evolution equations 
for the mode functions ( the function multiplying the cre- 
ation/annihilation operator) U^^ and V^^ ( i,j runs over 
1 and 2 for two field case and r for helicity ±) and cal- 
culate the occupation number of gravitino created from 
vacuum through these mode functions. In general, one is 
interested in the physical mass eigenstates {ipi,tp2) which 
are non-trivial c;onibinations ( with matrix coefficients) 
of gravitino 9 and matter chiral fermions T. So, strictly 
speaking, one would need to diagonalize the Hamiltonian 
at each moment of field evolution to keep track of the 
mass eigenstates and their abundance. This diagonaliza- 
tion process is rather involved^. Here we use a further 
simplification for our numerical analysis because wc are 
only interested in the asymptotic value of these abun- 
dances. Indeed, since the mixing is small at such late 
times, wc can simply follow the fields of interest {9,T). 
That is, at late times, these fields are approximate mass 
eigenstates and, therefore, their occupation numbers cor- 
respond to those of (x>'i, x>'2)- The validity of this approx- 
imation will be confirmed if we find that our occupation 
numbers cease to evolve at the time scales of interest. 

With this simplification, we may represent the mode 
decomposition in the following familiar form, 



X' (x) = 



(27r) 



3/2 



u^^k, 7?)<(fc) + y;^ (fc, v)bY{-k) 



(3.28) 



We then define the spinor matrix U- and U-^ 



^We refer the readers to [7] for the general discussion. 



v2 vy 
71 ^'-'^^^^ 



(3.29) 

with eigenvectors of the helicity operator cr • v/|v|, = 
^ Q j and ijj- = ^ " j . Using these spinor matrices, the 

field equation of motion (3.24) has a following simple 
form in terms of the matrices U+ and U- , 



a{t)U± 



(3.30) 



We can then expand U± in terms of positive and neg- 
ative frequency solutions, 



U+{t) ^ 1 



M 



MV'^ -if a,dt' 



1/2 



U-{t) 



1 



e 

1/2 



M 

LO 



1/2 



M 
LO 



1/2 



(3.31) 



where diagonal matrix lo = ^/k"^ + M^. a and (3 are pre- 
cisely the generalization of Bogolubov coefficients. In- 
deed, in the same way as Bogolubov coefficients, we can 
calculate the occupation number of i*'' fermion eigen- 
states from fi as 

Ni (t) = {f3*^f)^^ ( no summation for i ). (3.32) 

Wc also keep in our mind that, because of nontrivial mix- 
ing of fermion mass eigenstates for the case of coupled 
field system, we need an extra care about the identifica- 
tion of inflatinos and gravitinos. 

We solved the coupled mode equations (3.30) numeri- 
cally to obtain the occupation numbers for Ng and N^j. 
These arc plotted in Figure (1) as a function of comov- 
ing momentum at time 1000 in units of infiaton mass 
m<j ~ ^ which gives a typical time scale for infiaton 
oscillations. We have used the typical parameter values 
jl = 0.0001, A = 0.001 with an initial infiaton amplitude 
0.2Mp. The O'Raifeartaigh field $i has an initial ampli- 
tude /t^Mp and we normalized a{t) to be one at the start 
of our calculation. 
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FIG. 1. Gravitino abundance as a function of comoving 
momentum in units of inflaton mass. 



Ng in our plot corresponds to the field abundance 
whose mass converges asymptotically to gravitino mass 
\m\/Mp which shows up in the first element of the mass 
matrix M in cqn(3.26), and it gives an estimate of grav- 
itino abundances produced in preheating era. We point 
out that, in the presence of time dependent background 
as usually the case in dealing with cosmological problems, 
kinetic term of scalar fields can cause the supersymme- 
try breaking as we can easily see from the supersymme- 
try transformation of chiral fermion, /^^ ( superpartner of 
scalar x, say), 



7° dx 
V2 dt 



(3.33) 



where m-^ is defined in eqn(2.5). Therefore we should 
be aware that the values for A'^i and A^2 at intermediate 
time ( i.e. the time when a and are still far from its 
settlement in the vacuum) does not represent either of 
the inflatino or gravitino abundance because there exist 
non-negligible contributions of supersymmetry breaking 
from both inflaton and hidden sectors. 

The time when this SUSY breaking contribution of a 
and rria becomes comparable with that of supersymme- 
try breaking sector is beyond the range of our numerical 
integration. We checked, however, that in such a small 
parameter range with so small initial amplitude for as 
in our model, A^i and N2 converge to asymptotic values 
and do not change anymore at relatively early stage even 
when a still keeps its oscillation. This verifies that Fig.l 
should represent a good asymptotic behavior for grav- 
itino abundance. Especially, the cut-off scale [k ~ m^) 
of the comoving momentum for the number density does 
not change anymore. This reflects the fact that the oscil- 
lation scale of high momentum mode {k ^ rria) is much 
larger than that of the background field (m3/2), so high 



momentum modes behave adiabatically and do not re- 
sult in the non-adiabatic amplification anymore in pre- 
heating era. 

We further comment on the subtle problems in iden- 
tification of gravitinos and inflatinos in terms of mass 
eigenstates. Whatever field which dominates the local 
supersymmetry breaking is considered to be longitudi- 
nal components of gravitinos via super Higgs mechanism. 
The non-thermal production of fermionic fields is efficient 
just after inflation when the kinetic term of inflaton field 
governs the supersymmetry breaking in our model. This 
is nothing but the non-adiabatic field amplifications of 
longitudinal components of gravitinos which 'eats' the 
fermionic partner of inflaton field, i.e. inflMino. Hence 
when the fermion preheating is robust just after infla- 
tion, it is to-be inflatino, not to-be gravitino of our inter- 
ests, that is amplified via parametric resonance effects as 
longitudinal components of gravitinos free from Planck 
suppression. 

We also should mention that, because of the couplings, 
hidden sector may not be the sole cause for supersymme- 
try breaking even in the vacuum, and in fact, there could 
be still supersymmetry breaking from inflaton sector in 
the vacuum at a later time. It, however, can be shown 
that in the vacuum the supersymmetry contribution from 
inflaton sector is at most of order compared with 
\F\^ (~ jjL^) due to hidden sector [23,24], and we stiU ob- 
serve the dominant contribution of local supersymmetry 
breaking from the O'Raifeartaigh field at a later time. 
We thus expect our plot of Ng still represents a fairly 
good overall behavior of gravitino abundance in asymp- 
totic regime. 

For the comparison with the gravitino number density 
constraints from photo-dissociation process in nucleosyn- 
thesis for the case of thermal production of gravitinos in 
thermal bath (3.6), we need to integrate Ng{k) over the 
comoving momentum space. As usually the case with the 
preheating of fermions, our plot also indicates that occu- 
pation number as a function of comoving momentum k 
can be as large as of order unity at most up to the or- 
der of inflaton mass scale, kcutoff — ™ct and decreases 
exponentially for bigger k. So the number density for 
longitudinal components. 



^^3/2 



i^k^dk 



(3.34) 



during the preheating is at most 



A6 

n3/2(Ve) ^ kl^toff -ml- W'^'M/. (3.35) 
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We can now estimate the upper bound of the ratio of 
gravitino number density 713/2 to entropy density in anal- 
ogy with ( 3.4 ). 



^3/2 



S{tr) 



0.445,T| 



RH 



RH 



(3.36) 



and substitution of 1.66 • T^h\/9*/-^p H t^^ ior tr 
gives us the estimate of n^jils after reheating, 



^3/2 ^ n3/2(Ve)^pre^flH(l-66)^ 

s ~ 0.44M2 



(3.37) 



We expect this efficient gravitino production occurs well 
within the time range of typical oscillation of super- 
symmetry breaking field and we can substitute t^re ~ 
~ (a„/x)~^ ~ 10^M~^ and eqn (3.35) in above 
equation to obtain the upper bound, 

This upper bound of n/s for the gravitinos from non- 
thermal production in our model is thus smaller than 
eqn(3.6) of thermal scattering by at least two orders of 
magnitude. 

We therefore find that our model does not lead to the 
overproduction of gravitinos due to nonthermal process 
in preheating period, and reheating temperature con- 
straint due to this effect is less severe than that of grav- 
itinos produced by the scattering in thermal bath dur- 
ing the reheating period. We also point out that the 
expression given by eqn(3.37) was derived in a general 
setting and it can be used to obtain, in combination with 
eqn(3.6), the estimate for the relative significance of the 
gravitino production in thermal and non-thermal pro- 
cesses once the model of supergravity infiation is given. 



IV. CONCLUSION AND DISCUSSION 

We showed in this letter a realistic supergravity infla- 
tion model which breaks local supersymmetry in the vac- 
uum dominantly via f -term coming from O'Raifeartaigh 
field in the hidden sector. We emphasized the significance 
of radiative corrc;ction in supersymmetry breaking sector 
to evade the moduli problem, and subsequently obtained 
the non-minimal Kahler potential arising from this loop 
correction. Using this non-minimal form of Kahler poten- 
tial, we analyzed the possible non-thermal production of 



gravitinos in preheating era. The emphasis in our anal- 
ysis was on the longitudinal components of gravitinos 
which do not suffer from Planck suppression and hence 
potentially could lead to robust amplification via para- 
metric resonance effects. Wc showed that the comoving 
number density of fermionic mass eigenstate converges to 
its asymptotic value at an early stage of preheating era 
which is well before the time when the supersymmetry 
contribution comes from hidden sector fields. Physically, 
this indicates that the large comoving modes (fc > rria-) is 
much larger than the typical coherent oscillation of back- 
ground fields at later times {k ~ TO3/2 )• Hence comoving 
number density for big modes behave adiabatically and 
we do not expect the parametric amplifications at later 
times. We also discussed the subtle problems in identifi- 
cation of gravitinos and inflatinos due to the non-trivial 
mixing of fermionic mass cigenstates. To-be inflatinos 
are longitudinal components of gravitinos just after in- 
flation, and its role is replaced by fermionic partner of 
O'Raifeartaigh fields at later times. 

We estimated the upper bound of number density of 
non-thermally produced gravitinos by integrating out its 
comoving occupation number over momentum space. Be- 
cause of the small mass scale of inflaton field, the typ- 
ical momentum scale of produced gravitinos and con- 
scqiicntly the momentum space over which occupation 
number is integrated out turns out to be small as well. 
This leads to the relatively small number density of grav- 
itinos and wc showed that it gives less significant con- 
straint than that of gravitinos which are produced by 
thermal scattering. 

We point out that the cases involving three and more 
superchiral fields are rather involved. We can basically 
follow the formalism discussed in section IIIB, but we 
need additional care in interpreting the fermion fields as 
a superposition of mass eigenstates because we cannot 
completely gauge away one of fermion fields via unitary 
gauge as we can do in the two field case [6,7]. The case 
including the gauge interaction terms and the model of 
other supersymmetry breaking mechanism besides hid- 
den sector supersymmetry breaking are also to be exam- 
ined. 
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